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ABSTRACT 



A method for reducing DC offset from a receiver signal. The 
method includes jointly (i.e., simultaneously) estimating 
such DC offset and channel impulse response, and reducing 
the DC offset in accordance with the estimated DC offset and 
the estimate of the channel impulse response. 
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FIGURE 4 



Store the received burst samples r(k) 
K= 1....N 



200 



Joint MMSE optimization with respect to DC offset (A) and 
channel coefficient vector (h) resulting in coupled system 
of equations that are solved simultaneously 

A 0 = E jr(/c)^h H "b(/c)| 
\= EjZ(kjb H (k)y E<(b~(k)[r(k)-A 



O_/202 



Forward computed DC offset ( A 0 ) and channel coefficient vector ( h 0 ) 
to data receiver 34 (Fig. 3) 
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FIGURE 5 



Store the received burst samples r (n) 
n=1 to N 



200 



Calculate initial DC offset from received samples 
as the average value over the burst 



Subtract the initial DC offset and establish timing 



202 



204 



Calculate error function based on midamble (known) bits 
or decoded data bits 

e{n) = r{n)-h (n) b(n) -A( n ) 



Update DC offset value 
A{n +1 ) = A(n) + ne(n ) 



,208 



Update channel coefficient vector value 
h{n+1) = h(n)+nb(n)e*(n) 



/-V^210 
212 



Forward computed DC offset A(n+1) and channel coefficient vector h(n+1) 
to data receiver at each sample instant n 



214 





Proceed to data receiver decoding 
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FIGURE 6 



Store the received burst samples r(k) 
k=1 ... N 



300 



Calculate initial DC offset from received 
as the average value over the burst 



I 



samples ^^y^QZ 



Subtract the initial DC offset and establish timing /^/304 



Update channel coefficient vector value 
h ("+1)=TjX b (A) [ r (k) -/» (n)] * 

M K=1 



306 



Update DC offset value (calculated over midamble length with 
known bits or a tentative decision of information bits) 

1 *** 

/4(n + 1)=— 2, [r(k)-h(n+1) H b(A)l 

M K*1 




Forward computed DC offset A(P) and channel coefficient vector h(P) 

to data receiver 34 
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FIGURE 10 
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METHOD FOR CORRECTING DC OFFSETS 
IN A RECEIVER 

BACKGROUND OF THE INVENTION 

This invention relates generally to radio frequency receiv- 
ers and more methods for reducing DC offsets in such 
receivers. 

As is known in the art, received radio frequency signals 
are converted to baseband using various receivers. With a 
homodyne receiver, the received radio frequency signal is 
mixed with the local oscillator whose frequency is equal to 
carrier frequency of the received radio frequency signal to 
translate the carrier frequency to DC and thereby provide 
"direct conversion" of the modulation on the received radio 
frequency signal to a modulation at DC. Hence, a homodyne 
receiver is sometimes referred to as a direct conversion 
receiver. 

While such direct conversion receivers offer the advan- 
tage of reduced cost, DC offset associated with such receiv- 
ers presents a major problem for receiver performance. More 
particularly, DC offset is produced from the homodyning. 
The level of the DC offset may be significantly larger than 
the desired signal, i.e., modulation, to be demodulated. Thus, 
DC offset compensation techniques are typically required. 
To insure flexibility for different operating conditions, DC 
offset compensation can be part of the digital baseband 
portion of the receiver, preferably a part of the digital signal 
processing (DSP) firmware. One application of direct con- 
version receivers is in mobile or cellular communication 
systems. In such systems, the radio channels received sig- 
nals also suffer from intersymbol interference (ISI) caused 
by multipath radio propagation and. transmitter and/or 
receiver filtering. 

More particularly, various DC offset compensation 
techniques, for both analog and digital domain, have been 
suggested. With regard to the Global System for Mobile 
Communications (GSM) environment, these suggested tech- 
niques are largely dependent on the sources of DC offset 
since they result in different levels of DC offset compared to 
the desired signal level. As discussed in an article entitled 
"Design considerations for direct-conversion receivers," by 
B. Rezavi, published in IEEE Transactions on circuits and 
systems — II: analog and digital signal processing, June 
1997, pp 428-435, two major mechanisms causing DC 
offset in direct, conversion receivers are: Local Oscillator 
(LO) isolation to the receiver's Low Noise Amplifier (LNA) 
and mixer inputs; and interference leakage to the LO (i.e., 
self mixing). The level of DC offset in this case is much 
larger than the level of the desired signal. 

Various techniques have been suggested to remove DC 
offset generated by these two major mechanism. Included in 
these techniques are: AC coupling; Offset cancellation using 
capacitor; Sample mean (average) estimation; Adaptive DC 
offset compensation for burst mode operation; and 
Re-integration. 

AC coupling, requires high-pass filter with corner fre- 
quency 0.1% of the data rate, which is less that 270 Hz for 
GSM. Problems associated with this approach are elimina- 
tion of signal content around DC, group delay characteristic 
of the filter and settling time in Time Division Multiple 
Access (TDMA) environment. 

With offset cancellation using a capacitor for TDMA 
systems, the offset in the receive path can be stored on a 
capacitor during the idle mode and subtracted from the 
signal during actual reception in the burst. The major issues 
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with this technique are kT/C noise and problems when 
interfere r is stored along with other offsets. 

With the sample mean (average) estimation based 
technique, such technique includes averaging over sufiB- 

s ciently long period and subtraction, as described by A. 
Bate man and D. Haines, in "Direct conversion transceiver 
design for compact low-cost portable mobile radio 
terminals," in Proc. VTC'89, pp 57-62. In TDMA systems, 
averaging is usually performed over the burst duration. The 

10 issues associated with this approach are: it does not address 
DC offset changes within burst; and it may introduce some 
bias since burst does not have zero DC component due to 
different number of zeros and ones in the data stream. 
Besides its simplicity this method has some desirable sta- 

35 tistical properties. Sample mean estimate is optimal DC 
offset estimate in zero-mean noise with Gaussian probability 
density function in minimum mean-squared error sense (also 
minimum variance unbiased estimate) and maximum like- 
lihood sense. Even when the probability density function 

20 (pdf) of the noise is not known, the signal average is the best 
linear unbiased estimate, see Fundamentals of statistical 
signal processing: estimation theory, Prentice Hall, 1993 by 
S. Kay. 

Adaptive DC offset compensation for burst mode opera- 

25 tion is presented in a paper entitled "Adaptive DC offset 
compensation algorithm for burst mode operated direct 
conversion receivers," by S. Sampei and K. Feher, Proc. 
VTC92, pp 93-96. This approach utilizes the known bits 
from the preamble to acquire DC offset, typically 3-5 bits. 

30 Simulations have shown that the technique is efficient in the 
cases where the DC offset ratio (amplitude ratio of DC offset 
to the maximum amplitude of the transmitted symbol) is less 
than 40%. Degradation of the performance is up to 1.5 dB. 
A similar approach has been presented by J. Bergmans, 

35 Digital baseband transmission and recording, Kluwer Aca- 
demic Publishers, Section 8.8.2. in a more general commu- 
nication scenario of continuous reception. In essence it is a 
form of dynamic loop tracking DC offset. 

40 Digital compensation based on Least Mean Square (LMS) 
adaptive algorithm has been presented in by J. Cavers and 
M. Liao, in "Adaptive compensation for imbalance and 
offset losses in direct conversion receivers," IEEE Transac- 
tions on Vehicular Technology, November 93, pp 581-588. 

45 Such paper describes models and theoretical development of 
receiver and transmitter compensation (modulator and 
demodulator). The least mean square (LMS) algorithm has 
been applied to set the parameters of the compensating 
circuit (compensates for In-Phase (I) and Quadrature (Q) 

50 gain and phase imbalance and DC offset). Effectively the 
system has three adaptive coefficients. A drawback of the 
algorithm is the long convergence time and sensitivity to the 
selection of LMS step-size parameters. 

The re-integration approach is presented in a paper by B. 

55 Lindquist, M. Isberg and P.Dent, entitled "A new approach 
to eliminate the DC offset in a TDMA direct conversion 
receiver," In Proc. VTC'93, pp 754-757[7], The idea is to 
differentiate signal, digitize it and then re-integrate, thus 
eliminating DC component. It is based on adaptive delta 

60 modulation and since there is no time constant it targets 
TDMA direct conversion receivers. Simulation results pre- 
sented in the paper indicate signal to noise ratio (SNR) 
degradation of 1 dB in static channel for the Bit Error Rate 
(BER) range 1% to 0.1%. 

65 In a GSM system, the whole burst is stored and the 
all-digital techniques described above may be adapted to 
extract DC offset. Thus, referring to FIG. 1, the data receiver 
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stores the burst of data, r(k), where k-1 ... N and N is the 

number of samples in the burst. Each burst includes a e(k) - rik) - a - ^\^ n )' h (k - n)] 

mid -amble having a known sequence of bits disposed #> 
between data, (i.e., information bits) as shown. Such known 

sequence of bits is used to aid in equalization and more 5 ^ accordaQce ^ anomer fcature of the invention> a 

particularly for enabling computation of the channel impulse communication system is provided wherein information is 

response (CIR). As shown in FIG. 1, an estimate of the DC transmitted through a channel as a series of bursts, each 

oflset, A, is calculated. The estimated DC offset, A, where burst having a predetermined series of bits and a series of 

k«l ... N, is subtracted from the received burst. The result, information bits. The system includes a receiver for receiv- 

r(k)-A, where k-1 ... N, is processed to find an estimate of 10 ing the transmitted information. The receiver has a processor 

the channel impulse response (CIR), h. The estimate of the programmed to simultaneously solve the following equa- 

channel impulse response (CIR), h, can be obtained by tions from: (a) the predetermined series of bits; (b) a 

cross-correlating [r(k)-A] with the known mid-amble bit tentative decision of the information bits; or, a combination 

sequence. °f tne predetermined series of bits and the tentative decision 

of the information bits: 

The complexity of the DC oflset cancellation in GSM 

system is related to other signal processing functions per- $ f{e{k)) 

formed in baseband (synchronization, equalization). ~~ ^ =0 

Residual DC oflset may affect the performance of data „ „ , 

J K 20 df(e{k)) 
receiver. : — = 0 



dh 



SUMMARY OF THE INVENTION 



where: 

In accordance with the present invention, a method is 25 f(c(k)) is a function, usually quadratic, of the estimation 
provided for reducing DC offset from a received signal. The error, e(k), where 

method includes jointly estimating such DC offset and a e(k) is the difference between the received signal and an 

channel impulse response, h, and reducing the DC offset in estimate of the received signal; i.e., 

accordance with the estimated DC offset and the estimated 

channel impulse response, h. 30 c(Jt) , KA) _ A - £ [b[n)h(k-n)} 

In accordance with another feature of the invention, a " 
communication system is provided wherein information is 
transmitted through a channel having a discrete channel 

impulse response h(k), where k is a time index, to produce BRIEF DESCRIPTION OF THE DRAWINGS 

at an output of the channel, a signal, r(k), where: 35 These and other features of the invention will become 

more readily apparent from the following detailed descrip- 
r(*) = /4 + Y b{n)h(k-n) + N(k) l * on wnen taken together with the following drawings, in 

„ which: 

FIG. 1 is a block diagram of a data receiver in accordance 
40 with the PRIOR ART; 

where: FIG. 2 is a block diagram of a radio communication 

A is DC offset; system according to the invention; 

FIG. 3 is a block diagram of a data receiver used in the 
Y b(n)h(k-n) 45 rac ^° 00111 mumcat i° n system of FIG. 2 according to the 

^ invention; 

FIG. 4 is a flowchart showing the process used to jointly 
estimate DC offset and the channel impulse response for the 
is a modulated signal transmitted over the channel 42; system of FIG. 2 according to the invention; 
b(n) are transmitted data symbols; and 50 FIG. 5 is a flowchart showing the process used to jointly 

N(k) is additive noise. estimate DC offset and the channel impulse response for the 

The system includes a receiver for receiving the transmitted svstem °f FIG - 2 based ° n U ™ { Mean ^ uare error > san2 P le 
information. The receiver has a processor programmed to b y sam f^ over a bursl duratl0a according to one embodi- 
solve the following equations simultaneously: ment of tne invention; 

55 FIG. 6 is a flowchart showing the process used to jointly 
df^ic)) estimate DC offset and the channel impulse response for the 

^ " = 0 system of FIG. 2 based on Least Square Error, by processing 

over an entire burst in a specified number of iterations, 
d fM k )) _ o according to another embodiment of the invention; 

dh 60 FIG. 7 shows the probability density function (pdf) esti- 

mate of the DC offset estimation error for high signal-to- 
wnere . noise ratio (SNR), i.e., with no noise, according to the 

PRIOR ART (i.e., "sample mean") and according to the 
f(e(k)) is a function, usually quadratic, of the estimation invention (i.e., "joint estimate") performed over the entire 

error, e(k), where 65 burst; 

e(k) is the difference between the received signal and an FIG. 8 shows the probability density function (pdf) esti- 
estimate of the received signal; i.e., mate of the DC offset estimation error for a SNR=8 db, for 
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both the PRIOR ART method and in accordance with the 
invention over the entire burst; 

FIG. 9 shows the probability density function (pdf) esti- 
mate of the DC offset estimation error for a SNR=8 db, for 
both the PRIOR ART method and in accordance with the 5 
invention over the mid-amble portion of the burst; and 

FIG. 10 shows bit error rate (BER) as a function of Energy 
per bit/Noise density, (Eb/No), for no DC offset 
compensation, a compensation according to the PRIOR ART 
method, and in accordance with the invention over both the 
whole burst N (JE-N) and over the mid-amble portion of the 
burst M (JE-M). 



DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 



where k is a time index. 



10 



understand the process used to remove the DC offset, A, let 
it first be assumed that dynamic range of the converter is 
large enough to accommodate both DC offset and the desired 
signal (i.e., the modulation at baseband), and that the 
sampled version of the complex baseband signal at the 
output of the analog to digital converter (ADC) 34 can be 
represented as: 

r(k) = A + Y J - n) + N(k) 



where: 

A is the DC offset; 



15 



Referring now to FIG. 2, a block diagram of a commu- 
nication system 10, here a TD MA communication system is 
shown to include a transmitter, here represented by a trans- 
mitter 12 for producing a radio frequency carrier modulated 2Q 
by bursts of data bits. Each burst includes a known sequence 
of bits and a series of unknown, here information, bits. Such 
known sequence of bits is used to aid in equalization and 
more particularly for enabling computation of the channel 
impulse response (CIR). For example, with a GSM, each 25 
burst includes a mid-amble having a known sequence of bits 
interposed between a pair of series of information, bits. 
Here, the mid-amble is used in a manner to be described to 
reduce DC offset from a receiver signal by jointly estimating 
such DC offset and the channel impulse response, and 3Q 
reducing the DC offset in accordance with the estimated DC 
offset and the estimated channel impulse response. 

More particularly, the information, after passing through 
a data block generator 14, interleaver 18, burst formatter 20, 
modulator and analog to digital converter (ADC) 22 is 35 
transmitted by a transmitter 24 through space (i.e., a radio 
propagation channel 26) to a receiver 28. More particularly, 
the receiver 28 is a direct conversion receiver having a 
homodyning receive section 30, an ADC 32, a data receiver 
34 (to be described in detail in FIG. 3), a burst demodulator 40 
36, a de-interleaver 38, and a decoder 40. It is noted that a 
channel 42 is thereby provided between the output of the 
burst formatter 20 and the output of ADC 32. Such channel 
42 may be characterized by a discrete impulse response h(k), 



45 



The receive section 30 includes an antenna 43 for receiv- 
ing the transmitted radio frequency signal having a carrier 
frequency ft,. The radio frequency signal carrier is modu- 
lated with desired data at the transmitter, in any conventional 
manner. Here, the receiver is used in a Time Domain 50 
Multiple Access (TDM A) system. The antenna 43 is fed to 
a mixer 44, as shown. Also fed to the mixer 44 is a local 
oscillator signal produced by local oscillator (LO) generator 
46. Here, the frequency of the signal produced by the local 
oscillator is also 4. The output of mixer 44 is passed through 55 
low pass filter 48 to remove higher-order harmonics from the 
mixing process. The resulting baseband signal, here a train 
of bursts, produced at the output of the low pass filter is fed 
to the analog to digital converter (ADC) 32 where the 
modulation is converted to corresponding digital data. The go 
data is processed in the data receiver 34 which includes a 
Digital Signal Processor (DSP), as indicated. 

Referring to FIG. 3, the data receiver 34 is shown to 
include a burst storage, or memory, 50, used by a DSP 52 to 
remove the DC offset from the burst, and to provide equal- 65 
ization which is assisted by a joint estimate of both the DC 
offset, A, and the channel impulse response, h. In order to 



is the modulated signal received by the ADC; 
k is a time index; 
b(n) are data symbols; 

h(k) is the impulse response of channel 42 (i.e., from 
transmitter baseband to receiver baseband) and 
includes the modulator 22, the transmitter 24, the radio 
propagation channel 26, the R. F. receiver 30 and ADC 
32; and 

N(k) is additive noise. 

This formulation is quite general since it captures com- 
plex (i.e., In-Phase (1) and Quadrature (Q)) signal. 

The received signal includes three terms: (1) the DC 
offset; (2) the signal passed through the channel 42, i.e., the 
modulated signal; and (3) thermal noise. Intuitively, if we 
are able to remove effect of modulated signal, we will reduce 
overall noise in the process of DC offset estimation. 

The problem can be formulated as joint estimation prob- 
lem of unknown DC offset and channel 42 coefficients which 
represent the estimate of the channel impulse response. Such 
channel 42 coefficients are required for the data receiver 34 
equalizer 53 (FIG. 3) provided by the DSP 52. That is, we 
want to minimize a function of the error e(k) with respect to 
the estimated DC offset, A, simultaneously with the mini- 
mizing the error e(k) with respect to the channel impulse 
response, h, where e(k) is a function of the actual received 
modulated signal, r(k), and the estimated modulated signal, 
f, where k is the Lime index. That is, the DC offset and 
channel impulse response are jointly estimate. In order to 
perform this joint minimization of f(e), were solve the 
following two equations simultaneously: 



dA 
df(e{k)) 
dh 



= 0 



e(*) = /tf)-A-£ {b{n)h{k~n)] 



In a GSM system, for example, we constrain the estimate 
of h to a channel length to 5 coefficients based on physical 
channel profiles. In that case, the received signal can be 
expressed in the vector form, where the symbols: *; T; and, 
H are for conjugate, transpose and Hermitian (conjugate and 
transpose) respectively: 

K*W+"h w b'(*>fJV(Jt) 
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where: within the burst has to be calculated. Because it is extracted 

-* represents a vector; and fr° m mid -amble portion, it should be performed after DC 

offset removal. 

"h -{ho h Y Since the statistical average is not available in practical 

" ' L 5 system due to finite length of received samples (GSM burst), 

is the channel vector of length (L+l) assumed to be fixed there are several approaches to implement the optimal 

within a burst, and the corresponding data vector is solution. All of them are based on certain approximation of 



b(k)-{b(k)...b(k-L)-] r 



the statistical averaging operator. 
10 General Solution 



The simultaneous solution results in: 



Here, f(e(k)) is the Minimum Mean Square Error func- ^ . , . „ _ , t 

tion. The optimization criterion is the Minimum Mean- . ^ Process for reducing DC offset from a direct conver- 
Squared Error (MMSE), where symbol E{ } denotes statis- ™ n re ff ceiver 1S Permed by jointly estimating such 
tical average DC offset and channel impulse response (i.e. coefficients 

representing channel 42 (FIG. 2)) and reducing the DC offset 
2 in accordance with the estimated DC offset and estimated 

mmE{\e(k)\ }, channel impulse response. More particularly, referring to 

FIG. 4, the DSP 52 is programmed by a memory in the DSP 
52 in accordance with the flowchart shown in FIG. 4, to 
perform the following method on the received signal, r(k). In 
e(jt)-K£)-^-iT"ir(A)-JV(A) step 100, the received burst samples, r(k) are stored in a 

memory 55 of the DSP 52. Next, in step 102, a joint, i.e., 
Differentiating with respect to the DC offset value and the simultaneous, Minimum Mean Square Error optimization 
estimated channel impulse response and equating the with respect t0 the DC offset, A, and channel impulse 
obtained gradients to zero the coupled system of equations _» 

is obtained. That is: 25 response (i.e., coefficient vector, h (k)), is performed result- 

ing in a coupled system of equations that are solved simul- 

«*wf 0 taneously: 

^ h A<f*E{r{k)-h 0 "h(k)} 
dE\e{}cf_ 30 

dA h 0 -F4b(k)b"(k)}- l E{b(k)Wk)-A 0 r} 

Next, in step 104, the computed DC offset, Aq, and the 
computed channel impulse response, h 0 , (i.e., coefficient 
A< r E{r{k)-~h 0 H ~b(k)} vector ("f^)) are forwarded to the data receiver 34, FIG. 3. 

_ _^ _ _ Thus, referring to FIG. 3, the stored r(k) has the DC offset 

h 0 =£{ b (*) b "(*)}- l £{ b (kMk)-A Q ]*} thereof, A, removed by subtracting the computed DC offset, 

lliese equations provide optimal solution in MMSE sense V *™ * e st £f d ?\ * » 

(It is important to note that using maximum likelihood 40 sl S nal ™ e 15 fed to e ^ all ^ r 53 

criterion we would obtain the same solution since the along with the computed channel impulse response, h 0 , to 

additive noise is assumed to have Gaussian distribution). aid in equalization and ISI reduction. 

The following should be noted: Provided that wc have a 

statistical average (which is not the case in practice and Adaptive Solution Based on LMS 

which will be addressed later), the solution is intuitive: DC 45 „ . t F , t 

v , , . Having the gradient of the solution, we can replace the 

offset is calculated by averaging over the received samples u i - . • i ,u . 

. , . & . , . . , statistical average by single point approximation. In that 

once the channel impact is removed and the estimated , ,.° < .r e . . f u 

, . v rr. . x , . case, an adaptive algorithm for pint estimation can be 

channel lm pulse response (i.e, coefficients) are obtained summarized as* 
after removing DC offset from received samples. Because *, . . 

the equations are cross-coupled, we have to start the iterative 50 1 • Calculate initial values: 
process, i.e. estimate one of the parameters first and then 
proceed with the computation having in mind that DC offset 1 ^ 

can be much larger than the desired signal. An initial ~ Nfc* 

estimate of the DC offset can be average over received 

samples neglecting the channel impact. This is also intuitive: 55 
The statistical average of the modulated signal, here, for and 

example, a Gaussian Minimum Shift Keying (GMSK) N being the number of samples required for averaging 

received signal, is zero. In addition we can approximate (up to burst length). 

desired signal passed through fading channel as Gaussian 2 . Subtract A(0) from received samples and extract timing 



3. For n-0, 1, 2 . . , 



distributed (channel coefficients of fading channel are com 
plex Gaussian in European Telecommunication Standard 

Institute (ETSI) models). In that case, the optimal DC offset — w 

estimate in the zero-mean noise with Gaussian probability 4«M«)-h (") b («)-*(") 

density function (pdf) is given by received signal average, A(n+i)-A(n)+ye(n) 
which is optimal in MMSE and in the maximum likelihood 65 

sense. Even when pdf of noise is not known, the signal T(«+1)=1T(/j)V^ 
average is best linear unbiased estimate. The timing estimate 
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4. Subtract the DC offset value from received samples for 2. Subtract A(0) from received samples and extract tim- 

further processing, where is the Least Mean Square ing; 

(LMS) step-size parameter determined for each specific 3 p or n .Q ( i, 2 . . . calculate iterative solution 
system being designed though simulation of the sys- 
tem. From the computational point of view, this solu- 5 , ^ 
tion is very staple. It requires only one inner product % + 1) = ± V jB» W ± Y - *<„)]' = 
in error calculation. A potential problem for the LMS {*» frf J M fci 
solution is the long convergence time, which can be 

much higher than mid-amble length. Therefore, we — V wiiWAi-Atoii 

have to rely on tentative decisions feedback from 3Q 
receiver equalizer 53 (FIG. 3). On the other hand, the 

adaptive solution could track slow changes in the DC _ 1 - i 

ofifeet and the channel response within the burst. A{n + l) ~ m L r k) ' h(n + [) m \ 
Thus, referring to FIG. 5, the DSP 52 is programmed by 
a memory stored in the DSP 52 in accordance with the 

flowchart shown in FIG. 5, to perform the following method The DC offset estimate, A(n+1) ( can be calculated over 

on the received signal, r(k): the mid-amble length or entire burst length, N; and 

In step 200, the received burst samples, r(n), where n goes 4. Subtract the final DC oflset value from received 

from 1 to N, are stored in a memory provided in the DSP 52. samples for further processing. 

In the next step, 202, an initial DC offset, A(0), is calculated 2Q \ t [ s important to note that the index n in iterations does 

from the received samples as the average value over the no t correspond to time index as in LMS solution, but rather 

burst. In the next step 204, the initial DC offset, A(0), is it corresponds to the iteration step. Matrix inversion is 

subtracted from the received burst, r(n), where n-1 . . . N, avoided due to cross-correlation properties of the mid-amble 

(i.e„ r(n)-A(0). In the next step 206, an error function, i n the GSM system, so that the channel estimate is a simple 

e(n)-r(n)-T // (n)'b (n)-A(n) is calculated. In the next step 25 cross-correlation with known bit sequence. The number of 

208, an update of the DC offset value, A(n+1)-A(n)4^e(n) iterations required has to be determined via simulation, 

is performed. In the next step 210, an update the channel referring to FIG. 6, the DSP 52 is programmed by 

a memory stored in the DSP 52 in accordance with the 

impulse response (i.e., coefficient vector value) h (n+l)= flowchart shown in FIG. 6, to perform the following method 

"h (n)+jutT (n)e*(n) is performed. In the next step 212, the 30 0D the rece i ved signal, r(k): 

computed DC offset A(n+1) and channel impulse response 1° step 300, the received burst samples, r(k), are stored in 

(i.e., coefficient vector) h(n+l) are forwarded to the data * memory provided in the DSP 52. In the next step, 302, an 

receiver 34 at each sample instant n. The process returns to initial DC offset, A(0), is calculated from the received 

step 206 until the end of the burst, step 214. When the burst samples as the average value over the burst. In the next step 

is ended, the DSP 52 proceeds to data receiver 34 and 35 304, the initial DC offset, A(0), is subtracted from the 

decoding (i.e., decoder 40, FIG. 2), step 216. received burst. In the next step 306, the channel impulse 

Another possibility for the adaptive solution is to employ response (coefficient vector value) 
Recursive Least-Squares (RLS), however it is more com- 
putationally intensive. _ 1 v " 

Least-Squares Solution 4U M t-i 
The solution using the deterministic criterion of least- 
squares has the same form as the optimal MMSE solution are up d a ted. In the next step 308, the DC offset value 
except that the statistical average is replaced by a time ( ca i cu lated over the mid-amble length with known bits or 
average. Again, it is important to note that the DSP 52 is 45 lhe over the wnole burst with decocled bits> a tentative 
programmed to simultaneously solve the equations: decision of the information bits) 

dk ' A(n + 1) = — £[/W-Mn+l)"MA)] 

df{e{k)) = Q 50 * al 

8h 

are updated. In the next step 310, a determination is made as 

, , , , _ , , to whether the number of required iterations, P, been per- 

desenbed above from: (a) the predetermined series of bits; formed nave the Qumber of required iteralionS) P> been 

b) a tentative decision of the information bits; or, a com- 55 formed? If not> the rocess coa tinues to update the 

bination of the predetermined series of bits and the tentative estimated channel impulse response (Um coefficients) and 

decision of the information bits. 

The Least-squares soluti< 

1. Calculate initial value 



rn w , . , . , DC offset values, the process returns to step 306; if it has, the 

lae Least-squares solution can De summarized as: computed DC offset A(P) and channel impulse response 

(coefficient) vector h(P) is forwarded to the channel 
receiver, step 312. 



60 



1 f 



Simulation Results 



To quantify the performance improvement obtained by 
N being the number of samples required for averaging 65 performing joint estimation of DC offset and channel 
(up to burst length, N, to be determined by impulse response (i.e., coefficients), simulations have been 
simulation); carried out. The first set of simulations in Matlab compares 
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the probability density functions of the DC offset estimation 
error for sample mean estimator and joint estimator. The 
estimation error has been analyzed over a 1000 burst using 
a quasi-static channel modeling (fixed fading channel coef- 
ficients during the burst) and applying serial receiver analy- 
sis where DC offset has been added to the In-phase com- 
ponent of the signal after the signal de-rotation. Two 
scenarios have been analyzed: a high SNR (no noise) case 
and a low SNR=8 dB, and the results have been presented 
in FIGS. 7, 8, and 9 and TABLES I and II, below. 

TABLE I 

Statistical parameters of DC offset estimation 
error tor sample mean - averaging estimator (AE) and joint 
estimator (JE) performed over the whole burst N (JE - N) for 
high SNR fno noise') 





Static 


TU 


HT 


RA 


Mean - AE 


0.12 


0.001 


0.002 


0.0016 


Std - AE 


0.075 


0.1 


0.1 


0.1 


Mean - JE - N 


0.004 


le-4 


2e-4 


6c-5 


Std - JE - N 


0.045 


0.045 


0.05 


0.045 



10 



15 



20 



where Static refers to a direct connection in place of anten- 
nas for the radio propagation channel 26 (FIG. 2); RA is 
Rural Area; TU is Typical Urban Area; and, HT is Hilly 
Terrain, as defined by the ETSI. Further, "Mean" refers to 
statistical mean and Std refers to standard deviation. In high 
SNR case joint estimation eliminates the bias component of 
the DC offset residual in static channel, which is direct 
consequence of the modulated signal presence. At the same 
time in all channel conditions it reduced the variance of the 
estimate. The reduction in standard variation is about two 
times as quantified in Table I. 

In low SNR case standard deviation of both estimators is 
higher due to additional thermal noise. Joint estimator 
performed over the burst length eliminates bias for he static 
channel and reduces standard deviation by a factor of two as 
quantified in Table II. When joint estimation is performed 
over midamble length to reduce the computational complex- 
ity a few observations can be made. In all channel conditions 
standard deviation of the estimate has been larger since the 
sample size is smaller (midamble length versus whole 
burst). In static channel residual DC offset is eliminated, 
while variance is comparable to the sample mean case. In 
fading channels residual DC offset is two times smaller and 
standard deviation is 10% smaller when, using joint estima- 
tion. 



TABLE II 



Statistical parameters of DC offset estimation 
error for sample mean - averaging estimator (AE) and joint 
estimator (JE) performed over the whole burst N (JE - N) and 

the joint estimator performed over the mid-amble length 
M (JE - M). for SNR - 8 dB 





Static 


TU 


HT 


RA 


Mean - AE 


0.12 


0.004 


0.003 


0.004 


Std - AE 


0.08 


0.11 


0.11 


0.11 


Mean - JE - N 


0.005 


4e-4 


4e-4 


9e-5 


Std - JE - N 


.055 


.054 


.056 


.055 


Mean - JE - M 


.0001 


.002 


.001 


.001 


Std - JE - M 


0.097 


0.1 


0.1 


0.1 



25 



35 



45 
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1. The DC offset had random value up to 50 dB above the 
modulated signal with uniform distribution; 

2. The DC offset was constant during the burst duration; 
and 

3. The joint estimation technique has assumed the least 
complex implementation (i.e. single iteration for least- 
squares solution): 

(a) An initial DC estimate by averaging over the burst 

(b) A calculation of a 5-tap Channel Impulse Response 
(CIR) after initial DC subtraction 

(c) A modification of the DC estimate by subtracting 
the modulated signal is obtained by passing data bits 
from the entire burst (no errors in decision) through 
estimated CIR and averaging over the burst after- 
wards; and 

(d) A final Channel Impulse Response (CIR) estimation 
for the data receiver 34 is made after subtracting the 
DC estimate from the received signal. 

As shown in FIG. 10, the joint DC offset compensation 
performed over the entire burst length N (JE-N) has the 
potential of completely eliminating the performance degra- 
dation due to residual DC offset introduced by sample mean 
estimator. Performing the joint estimation over midamble 
length (JE-M) results in a loss of the order of 0.1 dB for SNR 
less than 8 dB, and no less for SNR larger than 8 dB. 

Other embodiments are within the spirit and scope of the 
appended claims. For example, while the invention has been 
described above for reducing DC offsets in a direct conver- 
sion receiver, the invention may be used with other types of 
receivers, such as heterodyne receivers. Further, the inven- 
tion may be used with systems other then GSM. 

What is claimed is: 

1. A communication system wherein information is trans- 
mitted through a channel having a discrete channel impulse 
response h(k), where k is a time index, to produce at an 
output of the channel, a signal, r(k), where: 

r(k) = A + Yj WM* - n) + N(k) 



where: 

A is DC offset; 



50 



55 



60 



is a modulated signal transmitted over the channel; 
b(n) are transmitted data symbols; and 
N(k) is additive noise, such system comprising: 

a receiver for receiving the transmitted information, 
such receiver having a processor programmed to 
solve the following equations simultaneously: 

df(c(k)) 



To quantify the impact on BER, Matlab simulations have 
been carried out for the static channel and summarized in 65 
FIG. 10. 

A few assumptions have been made, as follows: 



dA 
dh 



= 0 



where: 

f(e(k)) is a function of the estimation error quadratic, 
e(k), where 
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2. A communication system wherein information is trans- 
mitted through a channel as a series of bursts, each burst 
having a predetermined series of bits and a series of infor- 
mation bits, such system comprising: 

a receiver for receiving the transmitted information, such 
receiving having a processor programmed to simulta- 
neously solve the following equations from: (a) the 
predetermined series of bits; (b) a tentative decision of 
the information bits; or, a combination of the predeter- 
mined series of bits and the tentative decision of the 
information bits: 



df(e(k)) 

dA 
dfiejk)) 

dh 



where: 



= o 



10 



20 



f(e(k)) is a function of the estimation error, e(k), where 25 
e(k) is the difference between the received signal and an 
estimate of the received signal; 

e(k) = «k)-A-Y j [b{n)kk - «)]. 
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calculating an initial DC offset, A(0), from the stored 
received samples as the average value over the burst of 
samples; 

subtracting the initial DC offset, A(0), from the stored 
received burst to produce r(n)-A(0); 

calculating an error function, e(n)=r(n)- h "(n) 
"b(n)-A(n); 

updating the DC offset value, A(n+l)-A(n)+//e(n); 
updating a channel impulse response in accordance with 

lT(n+l)-.Tr(n)+//"b(n)e*(n); 
computing the computed DC offset A(n+1) and channel 

impulse response h (n+1) until the end of the burst. 
5. A method for determining DC offset in a receiver, such 
receiver receiving a burst comprising a known series of bits 
and information bits, such receiver comprising: 
storing the received burst samples, r(k), in a memory; 
calculating an initial DC offset, A(0), from the stored 

received samples as the average value over the burst; 
subtracting the initial DC offset, A(0), is subtracted from 

the stored received burst; 
updating a channel impulse response 



1 w 

A(n + 1) = [r{k)-A(n)Y 



3. A method for reducing DC offset from a receiver signal, 
comprising: 

storing received burst samples, r(k) in a memory; 35 
simultaneously computing: 



^o-£{K*)-h 0 "b(*)} 
'h 0 -£{"b(*)"b w (&)}- 1 £{"bCA;)[r(A)M 0 ]-} 



40 



to compute a DC of fset, Aq, and the channel impulse 
response, h 0 . 

4. A method for determining DC offset in a receiver, 45 
comprising: 

receiving burst of samples, r(n); 

storing the received burst samples, r(n) in a memory; 



updating a DC offset value over a from the known series 
of bits in the stored received samples, or over the entire 
burst, in accordance with: 



l M 



making a determination as to whether a required 
iterations, P, has been performed and, if not, the method 
continues to update the estimated channel impulse 
response and DC offset values; and if it has been 
performed, the computed DC offset and channel 
impulse response processed by the receiver. 
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